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Abstract 



In this small note we are concerned with the solution of Forward-Backward Stochastic 
Differential Equations (FBSDE) with drivers that grow quadratically in the control com- 
ponent (quadratic growth FBSDE or qgFBSDE). The main theorem is a comparison result 
that allows comparing componentwise the signs of the control processes of two different qgF- 
BSDE. As a byproduct one obtains conditions that allow establishing the positivity of the 
t> ■ control process. 
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1 Introduction 

This small note is concerned with forward-backward stochastic differential equations (BSDEs) 
in the Brownian framework, i.e. equations following, for some measurable functions b, a, f and 



<7, the dynamics 

Xl' x = x+ [ b(r,X*' x )dr+ [ a(r,Xp x )dW r , 
Jt Jt 

Y?> = g(X^) + f T f(r, X?,Y?*, Z^)ds - f Z, 

J s Jt 



where W a d- dimensional Brownian motion, (t,x) G [0, T] x R m and s € [t,T]. The function / 
is called generator or driver while g is named the terminal condition function. The solution of 
the FBSDE is the triple of adapted processes (X,Y,Z); Z is called the control process. 

In the last 30 years much attention has been given to this type of equations due to their 
importance in the fields of optimal control and finance. The standard theory of FBSDE is 
formulated under the canonical Lipschitz assumption (see for example |10j and references), but 
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in many financial problems drivers / which have quadratic growth in the control component 
appear i.e. when / satisfies a growth condition of the type \f(t,x,y,z)\ < C(l + \y\ + \z\ 2 ). The 
particular relation between FBSDE with drivers of quadratic growth in the control component 
(qgFBSDE) and the field of finance, stochastic control and parabolic PDE can be illustrated by 
the works [7J, [6], |10| and references therein. 

One of the fundamental results in BSDE or FBSDE theory is the so called comparison theorem 
that allows one to compare the Y components of of the solution of two BSDEs. In rough, given 
a terminal condition function g l , a driver f l and the corresponding FBSDE solution (X, Y l , Z l ) 
for i E {1,2}, if g 1 dominates g 2 and f 1 dominates f 2 in some sense then this order relation 
is expected to carry over to the Y components, i.e. Y 1 dominates Y 2 in some sense. Such a 
result is however not possible for the control components Z l . In this short note we give a type of 
comparison result for the control components Z, a so called comonotonicity result. This result 
allows one to compare the signs of the control processes Z 1 and Z 2 componentwise and as a 
side product one finds sufficient conditions to establish the positivity of the control process for a 
single FBSDE. 

This type of results can be useful in several situations, for instance in the numerics for such 
equations, since they allow to establish a priori heuristics that can improve the quality of the 
numerical approximation. This point of view is pertinent as the applications of FBSDE extend 
to the field of fluid mechanics (see [5]). 

A possible application of the results presented in this note lies in the problematic of showing 
the existence (and smoothness) of marginal laws of Y which are absolutely continuous with 
respect to the Lebesgue measure. This type of analysis involves showing the strict positivity 
of the Malliavin variance (in rough the Z component) of the solution of the FBSDE, (see e.g. 
[3]). The results in [3] were established for FBSDE whose driver function satisfies a standard 
Lipschitz condition in its spatial components and it is not possible to adapt the proof to cover 
the qgFBSDE setting of this work. 

From another point of view, the comonotonicity result is an interesting result in the con- 
text of economic models of equilibrium pricing when analyzed in the qgFBSDE framework. In 
such framework the equilibrium market price of risk can be characterized in terms of the con- 
trol process of the solution to a qgFBSDE. The difficulty is that the individual optimization 
problems underlying the characterization of the equilibrium requires the equilibrium volatility 
(the Z component of the solution to a certain qgFBSDE) to satisfy an exponential integrability 
condition as well as a positivity condition. Since the results of [3] cannot be applied or adapted 
to the qgFBSDE setting, the comonotonicity result presented here (and its corollary) provides 
conditions that ensure the positivity of the relevant process and hence may prove to be very 
useful in equilibrium analysis. An example of such type of problems can be found for example 
in 0. 

The results of this work originate in |4] where the authors give a comonotonicity result for 
FBSDE satisfying a standard Lipschitz condition and where the driver function is independent 
of the diffusion process X. In [14J the author extended the results of [1] to the qgFBSDE setting 
but was not able to include the dependence on X in the driver. The dependence of / in X is 
something that is quite common in the financial framework and that makes the applicability of 
|14| limited. This short note presents a full generalization of the results of |14j where the driver 
is now allowed to depend on X, this makes the conditions and analysis more involved but makes 
the result general enough that it can now be "broadly" applied to the standard financial setting 
where the driver / almost always depends on the underlying diffusion X. 

The note is organized as follows: In Section 2 we introduce some notation and recall some 
known results. The main results are then stated and proved in Section 3. 
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2 Preliminaries 

Throughout fix T > 0. We work on a canonical Wiener space (O, J-, P) carrying a d- dimensional 
Wiener process W = (W 1 ,--- ,W d ) restricted to the time interval [0, T] and we denote by 
T = (J 7 t)te[o,T] its natural filtration enlarged in the usual way by the P-zero sets. 

Let p > 2, then we denote by <S p (R m ) the space of all measurable processes (lt)j<=[o,T] with 
values in R m normed by H^Hsp = E[sup tG [ ,T] |^t| p ] 1 ^ p an d by 5°°(M m ) its subspace of bounded 
measurable processes. We also denote by T-L p (M. m ) the space of all progressively measurable 
processes (Zt)te[o,T] with values in R m normed by ||Z||%j> = E[( Jq \Z s \ 2 ds) p ^ 2 ] 1 ^ p . 

For vectors x = (x 1 ,--- ,x m ) G R m we write \x\ = (^2^1 (x 1 ) 2 ) 2 . V denotes the canonical 
gradient operator and for a function h(x, y) : R m x R d — > R we write V x h or V y h to refer to the 
first derivatives with relation to x and y respectively. 

We work with decoupled systems of forward and backward stochastic differential equations 
(FBSDE) for (t, x) G [0, T] x P m and s G [t, T] 

Xl' x = x+ [ b(r,X*' x )dr+ [ a(r,X*' x )dW r , (2.1) 
Jt Jt 

Y?° = giX^) + / T /(r,X^,y r ^,Z^)d S - C Z^dW r , (2.2) 

is Jt 

for some measurable functions b, a, g and /. 
We now state our assumptions. 

Assumption 2.1. The function b : [0,T] x PJ n -> K m and a : [0,7] x R m M mxd are 
continuously differentiable in space with derivatives uniformly bounded by a constant K and are 
^-Holder continuous in time, a is uniformly elliptic and |6(-,0)| and |<j(-,0)| are uniformly 
bounded. 

g : M m — >• R is bounded, continuously differentiable with bounded derivatives, f is a continu- 
ously differentiable function in space, uniformly continuous in the time variable and satisfies for 
some M > for all (t,x,y,z) G [0, T] x R m x R x R d , \ f(t,x,y,z)\ < M(l+ \y\ + \z\ 2 ) as well as 

\V x f{t,x,y,z)\ <M(l+|y| + |z| 2 ), \V y f(t,x,y,z)\ < M, \V z f(t,x,y,z)\ <M(l + \z\). 

Assumption 2.2. The spatial derivatives Vfe, Vex and Vg satisfy a standard Lipschitz condition 
in their spatial variables with Lipschitz constant K . 

Vy/ satisfies a standard Lipschitz condition with Lipschitz constant K and for all t G [0, T], 
x, x' G R"\ y, y' G R and z, z' G R d it holds that 

\V x f(t,x,y,z)-V x f(t,x',y',z')\ 

< K(l + \z\ + \z'\){(l + \z\ + \z'\)\x - x'\ + \y- y'\ + \z - z'\}, 
\V z f{t,x,y,z)-V z f(t,x',y',z')\ 

< K{{1 + \z\ + |z'|)|x - x'\ + \y- y'\ + \z - z'\}, 

The next theorem compiles several results found throughout [I], [8] and |14j . 

Theorem 2.3. Let Assumption \2. 1\ hold then for any p > 2 and (t, x) G [0, T] x R f/iere exists a 
unique solution Q f < x = (X l > x , Y f ' x , Z f < x ) of FBSDE -$ZM in the s P ace S p x S°° x W anc0 
f ZdW G BMO. 

1 BMO refers to the class of Bounded mean oscillation martingales, see [S] or for more details. 



3 



The variational process of Q t,x exists and satisfies for s G [t, T] 

V x Xl' x = I d + f V x b(r, X^VxX^dr + f V x a(r, X l ; x )V x X l ; x dW r , (2.3) 
Jt Jt 

VxY?* = V x g{X^ x )V x X^ + j\(Vf)(r, 8*'*), V^^da - ^ V x Z^ x dW r . (2.4) 

The triple Q t,x is Malliavin differentiable and its Malliavin derivatives are given by DQ t,x = 
(DX t ' x ,DY t ' x ,DZ t ' x ). The process {Z t s ,x ) s ^ t ^ has continuous paths, Z l ' x G S p and for <t< 
u < s <T the following representation holds 

D S Y*' X = Zj' x , F-a.s. and D U Y^ X = V x Y l s ' x {V ^^a^ X^ x ), F-a.s. (2.5) 

There exists a continuous function u : [0, T] x R m — > R such that for all (t, x) G [0, T] x M. m and 
s G [t,T] it holds that Y^ x = u(s,X l / x ) F-a.s.. 

Under Assumption \2.2\ the function u is continuously differentiable in its spatial variables and 
Z l f = (V x u)(s,X t s ' x )a(s,X t s ' x ) F-a.s. for all < t < s < T and x G M. m . 

Proof. Existence and uniqueness of the solution is quite standard either for the SDE (e.g. |13j ) 
or for the BSDE (see e.g. Theorem 1.2.12 and Lemma 1.2.13 in |14j). 

The variational differentiability and representation formulas as well as the path continuity of 
Z follow from Theorems 2.8, 2.9 and 5.2 in [8j (or Theorems 3.1.9, 3.2.4 and 4.3.2 of [13]). We 
emphasize that due to the continuity of the involved processes, the representation formulas (|2.5I) 
hold P-a.s. for all t G [0, T] and not just P (g> Leb-a.a. 

Lastly, the Markov property of the Y process is rather standard (see Theorem 4.1.1 of [14]). 
The differentiability assumptions on the driver and terminal condition function (Assumption 12. 2[) 
ensure that the function u is continuously differentiable in the spatial variables. A detailed proof 
of this can be found either in Theorem 7.7 in [2j or Theorem 4.1.2 in [14] , □ 



3 A comonotonicity result for quadratic FBSDE 

In this section we work with a ci-dimensional Brownian motion W on the time interval [0, T] for 
some positive finite T. Throughout let (t, x) G [0,T] x M m . Our standing assumption for this 
section is as follows. 

Assumption 3.1. Let Assumptions \2. il and \2.2\ hold. Assume that m = 1 and d > 1 

Remark 3.2. We note that it is possible to write the results of this section for multidimensional 
SDE systems (i.e. when m > 1) under the assumption that a is a square diagonal matrix and 
the system of forward equations is fully decoupled. There are many applications where such an 
assumption takes place (e.g. Jdj/). We write these result with m = 1 to simplify the presentation 
of this short note. 

For each i G {1,2} we define the SDE flU) with 6j and a { and BSDE ([22]) with terminal 
condition and driver given by gi and /j. We denote the respective solution of the system by 
^ x t,x,i^yt,x,i^ z t,x,ij^^ valued inlxMxl 1 ' for (t,x,i) G [0,T] x M x {1,2}. 

We define the vector-product operator, "0", as : M d x R d — > M. d such that 

a b = (ai&i, . . . , a d b d ), for any a = (oi, • • • , a d ), b = (6 l3 ■ ■ ■ , b d ) G R d . (3.1) 

With the convention that a b > means that for each i G {1, . . . , d}, Oj6j > 0. 

The aim of this section is to explore conditions such the following statement holds 

z t,x,i Q z t,x,2 y Q) p _ a ^ for any ^ x ) G x M and s G [t, T\. 
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Definition 3.3 (Comonotonic functions). We say that two measurable functions g, h : R — y R 
are comonotonic if they are monotone and have the same type of monotonicity, i. e. if g is 
increasing or decreasing then h is also increasing or decreasing respectively. We say that g and 
h are strictly comonotonic if they are comonotonic and strictly monotonic. 

We now state our main theorem. 

Theorem 3.4. Let Assumption^ hold and for (t,x) G [0,T] x R define (X t ' x ' i ,Y t > x > i , Z^) 
as the unique solution of FBSDE \2. l\) - [2Jfy for i G {1,2}. Suppose that x i— > g%{x) and x i— > 
fi(-,x,-,-) are comonotonic for alii G {1,2} and further, that gi,g 2 are also comonotoni^. If it 
holds for all s G [t,T] that 

a 1 (s,X t s ' x ' 1 )Qa 2 (s,X t s ' x ' 2 )>0, P-a.s., (3.2) 

then 

z t,x,l z t,x,2 > Q) P . 0iSij f or any G [ ,T] x R and s G [£,T]. (3.3) 

^/Sl, 52 are strictly comonotonic and inequality (|3 . 2[) /toWs strictly then (|3 . 3|) is a/so strict. 

Proof. Throughout take £ G [0,T], x G R and let i G {1,2}. According to Theorem 12.31 
for each i G {1,2} there exits a measurable deterministic, continuously differentiable func- 
tion (in its spatial variables) m : [0, T] x R — > R such that Yg' x ' 1 = Ui(s, Xt' x ' 1 ) and zt' x ' 1 = 
(V x Ui)(s, X t s ' x,l )a{s, Xt' x ' 1 ) P-a.s. We have then P-a.s. that for any s G [£, T] (recall that <jj is a 
vector and Vitj a scalar) 

z t,*,l Q z t,*,2 = ( ( v xUl )( S( ^,l) a^A^' 1 )) ((V^ 2 )(s,A^ 2 ) a 2 (s,A^ 2 )) 

= ((Ti(s, A*' 2 '' 1 ) (T2(s, A*'^' 2 )^ (V x ui)(s,X t s ' x ' 1 )(V x u)(s, X t s ,x,2 ). (3.4) 

A standard comparison theorem for SDEs (see [13J) yields that for any fixed £ and T the mappings 
x h- > A^'* are increasing. This, along with the fact that g± and g 2 are comonotonic functions, 
implies that for fixed £ and T it holds that x i— > gi(Ay X>1 ) and x i— > ^(A^' 2 ) are a.s. comono- 
tonic. A similar argument implies the same conclusion for the drivers /j, i.e. a; h- > A*' 21 ' 1 , •, •) 
and x H> /2(-, A.*' a ' 2 , •, •) are a.s. comonotonic. 

Using the comparison theorem for quadratic BSDE (see e.g. Theorem 2.6 in |12j ) and the 
monotonicity (and comonotonicity) of x i— > <7j(A~ 2 '*) and x h- > A.*' 2 '*, ■, •) we can conclude 
that x i— > F 1 ' 1 ' 1 is also a.s. monotone. Furthermore, since x H )• (^(A™ 1 ' ), x i— >■ ^(A^' 2 ), 
x i—T- /i(-, X t,x, ' L -, •, •) and x h- > /2(-, A.*'^' 2 , •, •) are comonotonic the same comparison theorem 
yields that the mappings x i— > Y 1 ' 01 ' 1 and x i— > y*' a ' 2 are also a.s. comonotonic. Equivalently, one 
can write for any (£,x) G [0,T] x R that (notice that Vu exists according to Theorem I2.3P 

((V«wi)(t, x), (V^ 2 )(£, x)> > 0. (3.5) 

Therefore, combining (|3.5p with (|3.2p in (j3.4j) we easily obtain 

Zl' x ' l {uj) Z^' 2 ^) > 0, P-a.s. w G fi, (£, x) G [0, T] x R, s G [£, T]. 

Under the assumption that g\ and 52 are strictly comonotonic it is clear that inequality (|3,5p is 
also strict. Furthermore, if one also assumes that the inequality in (|3,2p holds strictly for any 
(£, x) G [0, T] x R then fl3J3]) also holds strictly. □ 

2 This implies that x \-¥ /].(•, x, ■, •) and a; h-s> / 2 (-, a;, ■, ■) are comonotonic as well. 
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Unfortunately it doesn't seem possible to weaken the assumptions of the previous theorem. 
The key factor is the representation of Z l ' x via the function Y*' x = u(t, x) which needs to be 
continuously uniformly differentiable in the spatial variable and for that one needs Assumption 
2J to hold. 

We obtain an interesting conclusion of the previous result if we interpret the forward diffusion 
of the system as a backward equation. In terms of applications (as mentioned in the introduction) 
it is the next result that gives a condition that allows the user to conclude the positivity or 
negativity of the control process. 

In the next result we focus on just one FBSDE so we fix i = 1 and we omit this index. 



Corollary 3.5. Let the assumption of Theorem \3.4\ hold (fix i = 1). Take (t,x) £ [0,T] x R and 
let (X, Y, Z) be the unique solution of the FBSDE 



X t = x+ I b(s,X s )ds + I a(s,X s )dW s 
Jo Jo 



(3.6) 



Y t = g(X T ) + I' f(s, X S ,Y S , Z s )ds - F Z s dW s . (3.7) 
Jt Jt 

Then, if x i— > g(x) and x i— >• f(-,x, ,-,-) are increasing (respectively decreasing) functions, 
then Zt © o~(t,Xf) is F-a.s. positive (respectively negative) for all t € [0, T]. In particular, if 
the monotonicity of g and f (in x) is strict and if a is strictly positive then Z is either strictly 
positive or strictly negative (according to the monotonicity of g and f). 

Proof. Throughout let x € R and t £ [0, T\. We prove the statement for the case of g(x) and 
f(-,x, -j ■) being increasing functions (in the spatial variable x) and we give a sketch of the proof 
for the decreasing case. Rewriting SDE fj3.6|) BSDE leads to X t = X T - // b(s,X s )ds - 
jF a(s, X s )dW s . In fact we can still rewrite the above equation in a more familiar way, namely 

l-T 



Y t = g(X T )+ [ f(s,Y s )ds- ! Z s dW s , (3i 
Jt Jt 



it 

where Z s = a(s, X s ) for s £ [0, T], g{x) = x and f(t, x, y, z) = —b(t, y). 

At this stage we need to clarify the identification Z. = cr(-,X.). Let us write explicitly the 
dependence on the parameter x of the solution (X, Y,Z) of the FBSDE (|3.6p . (j3.8|) . i.e. we 
write (X, Y,Z) to denote (X, Y,Z). Note that the solution of the BSDE (|3,8p is the solution 
of SDE (|3.6p which is a Markov process. We can then write Y. = X. = u(-,X.) where it is the 
identity function (infinitely differentiable). Under Assumption 12 . 1 1 both X an Y are differentiable 
as a functions of x (see Theorem I2.3p . we have then Z. = (V x u)(-,X.)a(-, X.). And since u is 
the identity function with derivative being the constant function 1, it follows immediately that 
Z.=a(-,X.). 

Our aim is to use the previous theorem to imply this result. So we only have to check that its 
assumptions are verified. Comparing the terminal conditions of (|3.7p and (j3.8|) . i.e. comparing 
x i— > g(Xj,) with x i— >• g{Xj) = Xj, it is clear that both functions are almost surely increasing. 
Further, the driver function / of BSDE (|3.8p is given by f(t,x,y,z) = f(t,y) = —b(t,y) which 
is independent of x. Clearly x >— > f(-,x,-,-) and x i— > /(•, x, ■, ■) are comonotonic. 

Theorem 13.41 applies and we conclude immediately that 

Z t QZ t = Z t Qa(t,X t ) >0, P-a.s. te[0,T\. (3.9) 
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For the other case, when g is a decreasing function, the approach is very similar. We rewrite 
the SDE (|3.6p in the following way, 

-X t = -X T + f b(s,X s )ds- f [-a(s,X s )]dW t , t€[0,T\. 
Jt Jt 

The terminal condition of the above BSDE is given by x i— >• g(x) = —x evaluated at x = Xt 
and the driver f(t, x, y, z) = f(t, y) = b(t, —y) which is independent of x. Since g is a decreasing 
function, we obtain our result by comparing the above BSDE with (|3.7I) and applying the previous 
theorem. □ 

The above corollary allows one to conclude in particular the strict positivity of the control 
process. If one is only interested in establishing positivity (ignoring strictness), then one can 
indeed lower the strength of the assumptions. 

Lemma 3.6. Let Assumption \2. 1\ hold and m = 1. Assume that x h4 g{x) and x i— > f(-,x, 0, 0) 
are both monotone increasing then Zf a(t,Xt) > 0, P-a.s. for all t € [0, T], If x ^ g(x) and 
x i— > /(•, x, •, •) are both monotone decreasing then Z% © a(t, Xt) < 0, P-a.s. for all t € [0, T]. 

Remark 3.7. Again, as in Remark VS. "A it is possible to state and prove the same result for 
m > 1. One needs to impose that a is a square diagonal matrix and the SDE to be a decoupled 
system. 

Remark 3.8. It is possible to weaken the assumptions of this lemma as was done for Theorem 
4-3.6 in \1J$ or Corollary 2 in [9] '. Namely, the conditions are weakened to Lipschitz type con- 
ditions with the appropriate Lipschitz "constant", then one argues similarly but combining with a 
regularization argument. 

Proof. Throughout let t E [0, T] and x € R. Then due to the representation formulas in (|2.5p we 
have P-a.s. that 

Z t a(t, X t ) = D t Y t a(t, X t ) = V^V^)- 1 ^*, X t ) a(t, X t ). (3.10) 

It is trivial to verify that cr(t,Xt) o~(t,Xt) > 0. It remains to establish a result concerning the 
sign of V X Y and (V x X)~ l . 

Under the assumptions it is easy to verify that the solution of (I2.3|) is positive. The solution 
of V X X is essentially a positive geometric Brownian motion with a nonlinear drift and volatility 
which in turn implies that (V X X) _1 is also positive. If we manage to deduce a result concerning 
the sign of V X Y we are then able to obtain a weaker version of Corollary 13.51 

The methodology developed to deduce moment and a priori estimates for quadratic BSDE 
and illustrated in Lemma 3.1 and 3.2 of [8] (or Chapter 2 in [14J) allow the following equality 

V x Y t = E § [e T (e t )- l V x g(X t )WX T + j\e r et\V x f)(r, X r , 0, 0)V x X r ]dr\ F t ] , (3.11) 

where the process e and the measure P are defined as 

/ f* f(r,X r ,Y r ,Z r ) - f(r,X r ,0,Z r ) ^ X 
e t = exp| ly P? 4odr|, 

and P is a probability measure with Radon-Nikodym density given by 
dP M c ( /■ T /(r,X r ,0,Z r )-/(r,X r ,0,0) vt 
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Both (et)t£[o,T] an( i M are well defined. The first because y h- > f(-,-,y, •) is assumed to be 
uniformly Lipschitz and hence e is bounded from above and below and away from zero. The 
second follows from a combination of the growth assumptions on V«/ and the fact that J ZdW 
is a bounded mean oscillation martingalqfj (BMO). 

We have already seen that VX is positive and it also trivial to conclude that the process e 
also is. Given that g and / are differentiable, then saying that these functions are monotonic (in 
x) boils down to making a statement on the sign of (V x g){x) and (V x f)(-,x, 0, 0). If one assumes 
that g and _f(-,x,0,0) are monotone increasing in x then (Vg)(x) > and (V x /)(-, x, 0, 0) > 
for all x. Hence from (|3.1ip (and the remarks above) we conclude that V X Y is also positive. 
Returning to (|3.10p we have then that Zt cr(t, Xf) > which proves our result. 

The arguments are similar for the case when g{x) and f(-,x, 0, 0) are decreasing functions. □ 
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